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We present results of VLBI observations of the water masers associated with 
IRAS 4A and IRAS 4B in the NGC 1333 star-forming region taken in four epochs 
over a two month period. Both objects have been classified as extremely young 
sources and each source is known to be a multiple system. Using the Very Long 
Baseline Array, we detected 35 masers in Epoch I, 40 masers in Epoch II, 35 
in Epoch III, and 24 in Epoch IV. Only one identified source in each system 
associates with these masers. These data are used to calculate proper motions 
for the masers and trace the jet outflows within 100 AU of IRAS 4A2 and IRAS 
4BW. In IRAS 4A2 there are two groups of masers, one near the systemic cloud 
velocity and one red-shifted. They expand linearly away from each other at 
velocities of 53 km s _1 . In IRAS 4BW, masers are observed in two groups that 
are blue-shifted and red-shifted relative to the cloud velocity. They form complex 
linear structures with a thickness of 3 mas (1 AU at a distance of 320 pc) that 
expand linearly away from each other at velocities of 78 km s _1 . Neither of the 
jet outflows traced by the maser groups align with the larger scale outflows. We 
suggest the presence of unresolved companions to both IRAS 4A2 and 4BW. 

Subject headings: ISM:jets and outflows — masers — shock waves — stars:mass 
loss — stars :pre-main sequence 



1. Introduction 



Water masers are excellent probes of astrophysical flows. Their large flux densities and 
exceedingly compact sizes make them perfect Very Long Baseline Interferometer (VLBI) 
targets. Many masers observed in star-forming regions are thought to form in shocks either in 
or along the outflow of material commonly seen associated with young stellar objects (YSOs). 
They provide a high resolution probe of the base of the stellar wind that is unaffected by 
extinction from dust or extensive interaction with ambient material. The milliarcsecond 
resolution of VLBI is critical to resolve the individual maser components and outflows, 
particularly in high stellar density environments. 

Long term total power monitoring has shown that masers around low mass YSOs are 
more episodic than those around higher mass YSOs, however, their detectable phases are suf- 
ficiently long for pr oper motion measurements if the observations are spaced by no more than 



two t o three weeks (IFuruya et al.ll2003l ; iBrand et al.ll2003l ; IClaussen et al.lll996l ; IWilking et al. 



19941 ) . Furthermore, if they arise in the warm (400K) spatially confined po s t-shock region, 



as postulated in various mode ls (jElitzur. Hollenbach. fc McKed Il992l . Il989l ; iMacLow et al. 



1994 ; IMacLow &: Elitzurlll992l ). they should have space velocities sufficient to produce mea- 
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surable proper motions over one to three weeks. However, regular monitoring of low mass 
YSOs is necessary to ensure that VLBI observations are conducted during periods of high 
maser activity. 

Using the NRAO's Very Long Baseline Array (VLB A), proper motion studies have 
been made of water masers associated with high mass a nd intermediate mas s YSOs such 
as Cepheus A H W2 jTorrelles et all boOlh . S 1 06 FIR jFuruya et all lioooh. NGC 2071 
rtSeth et alJliooj ). and IRAS 20050+2720 MMS1 JFuruva et al.ll2005h as well as with lower 



mass YSOs such as Serpe ns SMM1, RNO 15-FIR, and IRAS 05413-0104 flMoscadelli et al. 



20061 ; IClaussen et al.lll998l ). For example, water masers in IRAS 05413-0104 (associated with 
HH 212 at a distance of 450 pc), were found to lie in a structure of 10 AU in length and less 
than 0.5 AU thick. The masers appeared to arise in shock-related structures which showed 
proper motions along the axis of the outflow of 60 km s^ 1 and displayed coherent structures 
over time scales of two to three weeks (but which was less discernible over time scales as long 
as several months). Comparison with molecular and infrared observations of IRAS 05413- 
0104 clearly demonstrated that proper motion studies of water masers in the region (within 
40 AU) of the central source are among the best tools for studying the kinematics of the jets 
emanating from embedded YSOs. 

IRAS 4, comprised of multiple sources whose submillimeter-dominated 'Class 0' spectral 
energy distributions suggest extreme youth, is located in the NGC 1333 star-forming region. 
As a nearby active and multiple young object, it presents a prime cand idate for VLB! 



obser vations of its associated water masers. Located at a distance of 320 pc fide Zeeuw et al. 
19991 ). the NGC 1333 molecular cloud hosts a double infrared cluster of about 200 YSOs 



ident i fied in near-infrared surveys (jStrom et al.lll976l; lAspin et al. 
1996; IWilking et al.l 120041 ) and x-ray surveys ]preibischlll997i : 



1994 



Lada, Alves. fc Lada 



Getman et al 



2002 



Preibisch 



20031 ). Far- infrared observations have not only revealed the higher luminosity sources in NGC 
1333, but also YSOs in the earliestphase of evolution such as the Class sources IRAS 4 

35 



and IRAS 2 flHarvev. Wilking. fc Jo\ 



1984 



Jennings et al.lll987l) . IRAS 4 was found to be a 



binary at submillimeter wavelengths (jSandell et al.lll99ll ); hence IRAS 4A and 4B. IRAS 4A 
is located 31" to the NW of 4B at a posit ion angle of -45° . Subs e quent ly, both IRAS 4A 
and 4B were found to be binary systems. iLay. Carlstrom. fc Hillsl ( 119951 ) performed single 
baseline interferometry in the submillimeter and found that 4A was a binary with an angular 
separation of 1.8" (or 580 AU at 320 pc). Interferometric observations in the millimeter 
and radio continuum co nfirmed the separation and determined that the posi tion angle of 
separation was « -20° flLoonev. Mundv. fc Welch! I2OO0I : iReipurth et~aD booj ). While the 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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easternmost source (IRAS 4A1) dominates the millimeter and radio continuum emission, the 
weste rnmost source IRAS 4A2) appears to be relatively more evolved with a warm ammonia 
core (IWootten fc Mangum Il993l : IShahl 120001 ) and a wealth of complex organic molecules 
(Bottinellii 2007, priv. comm.). S ubmillimeter and millime t er wa ve continuum observations 
have confirmed the prediction of [Lay. Carlstrom. fc Hilla (119951 ) that IRAS 4B is also a 
binary with components (4BE and 4BW) separated by about 10" (3200 AU in projection) 
along an E-W direction JSmith et al.lboool : lLoonev. Mundv. fc Welch! I2OO0I : ISandell k Knee 
2QQlh. § A s in the case of IRAS 4A2, complex organic molecules have been detected toward 
IRAS 4B (jBottineUi et al.l 120071 ). presumably from a hot core as sociated with 4BW. We 
note that IRAS 4BE has yet to be detected in the radio continuum (IReipurth et al.ll2002l ) or 
demonst rate any sign of outflow a ctivity. Infall motions have been detected toward IRAS 4A 
and 4B (IDi Francesco et al.ll200ll ). unambiguously toward the former, and strongly indicated 
toward the latter. 

The detection of water masers toward IRAS 4A and 4B suggests that they are YSOs also 
as sociated with mass ou tflow. Water masers we re first detected in sin gle dish observations 
by lHaschick et al.l (119801 ) and later monitored by lClaussen et al.l (119961 ). These observations 
established the maser emission as highly variable on monthly time scales, sometimes being 
completely absent while at other times reaching peaks of 10 Jy. Detected emission was found 
from —10 to +15 km s _1 but the angular resolution was not sufficient to determine which 
masers were associated with 4A and 4B. Beginning in 1983, VLA o bservations have shown 



water maser a c tivity associated on ly with IRAS 4A2 and/or 4BW ([Rodriguez et al.l 12002 



Woottenl Il998l ; iFuruya et al.l 120031 ). VLBA obse rvations of the water masers in IRAS 4A 
and 4B, acquired in 2003, have been presented by lDesmurs et al.l (120061 ). While masers were 
detected in only two of the four epochs observed, they confirmed the association of masers 
with IRAS 4A2 and expansion motion between the two spots was detected. In IRAS 4BW, 
the five maser spots detected were red-shifted and formed a chain 80 mas in extent with 
some masers displaying proper motion toward the north. 

Molecular outflows have been mapped toward both IRAS 4A and 4B. IRAS 4A is as- 
sociated with a highly collimated o utflow about 20,000 AU in extent seen in HCN a nd SiO 
with a position angle close to 20° iGirart Crutcher fc Raol Il999l : Ichoil I2OO1L l2005h . The 



origin of the outflow is likely IRAS 4A2 (jChoi 



20051 ). On a larger scale (« 1 x 10 5 AU) 



the outflow defined by CO and molecular hydrogen has a position angle of 45° (IBlake et al. 



19951 ; IChoi et al.l 120061 ) . The shift in position angle from small to large scale is perhaps due 



2 We adopt the naming convention introduced by 



Sandell fc 



^nee _j 200lh . The western (eastern) 



source is referred to as IRAS 4BI (IRAS 4BII) by ISmith et al 



Loonev. Mundv. fc Welch! |2000h . and as IRAS 4B (4B') bv |Di Francesco et al.ll200lh . 



(1200(1. as IRA S 4B (IRAS 4C) by 
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to a combination of a shift in the direction of t he magnetic field , an encounter with denser 
ambient gas, and precession of the outflow axis ( Choi et al. 20061) . A more compact outflow 



is associated with IRAS 4BW, with a position angle close to ° (IBlake et al 



with perhaps a second outflow with a position angle of -35° (IDi Francesco et al. 



1995 



Cho JboOlh 



200lh . 



In this paper we present VLBA observations of water masers associated with the IRAS 4 
region obtained over four epochs in 1998 and spaced by about one month. We also describe 
Very Large Array (VLA) observations, taken one month in advance of the start of the VLBA 
observations, that showed IRAS 4 to be a very active maser source and yielded absolute 
positions for the masers and the 1.3 cm continuum sources IRAS 4A1 and IRAS 4BW. With 
the VLBA, we detected 35 masers in Epoch I, 40 masers in Epoch II, 35 in Epoch III, and 
24 in Epoch IV associated with IRAS 4A2 and IRAS 4BW. We use these data to reveal 
the structure of shock fronts associated with stellar winds from these YSOs and to estimate 
proper motions for masers detected in all four epochs. The origin of the maser emission and 
its relationship to the stellar winds and larger scale molecular outflows from these YSOs is 
discussed. 



2. Observations and Reductions 

2.1. VLA Observations 

The VLA was used on August 12, 1998 to observe the water masers toward IRAS 4A 
and IRAS 4B to check that the masers were strong enough to observe subsequently with 
the VLBA. The array was in the B configuration. The phase tracking center for the array 
was set to be midway between the IRAS 4A and IRAS 4B positions, spread over a range of 
hour angles to optimize the u, v coverage. Two observation modes were used. In the first, 
observations were made with both circular polarizations, a total bandwidth of 1.5625 MHz 
(total velocity extent of 21 km s _1 ) and a correlator setup which gave 127 spectral channels 
each with width of 12.2 kHz (0.16 km s _1 ). In this mode the water masers were observed 
with reasonably high spectral resolution in order to spectrally resolve the masers, and the 
observations lasted for ~42 minutes (on source), 

In the second mode, the correlator was configured for continuum observations with one 
polarization pair (of IFs) set to a center frequency of 22264.9 MHz and a bandwidth of 50 
MHz, while the other polarization pair were set to a bandwidth of 0.781 kHz, and the center 
frequency adjusted so that the strongest masers were placed at the center of the narrow 
band. This setup allowed the phases to be tracked using the strong signal in the narrow 
band, in order to calibrate tropospheric phase changes to apply to the 50 MHz continuum 
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band (away from the maser emission). In this mode the continuum emission from the IRAS 4 
region was observed for ~3.5 hours (on source). 

The data were edited, calibrated, and imaged with the NRAO Astronomical Image 
Processing System (AIPS, Greisen 2003) in the standard manner. In the high spectral 
resolution data, Stokes I images were made using uniform weighting of the u,v data, of each 
spectral channel, and were CLEANed to an rms noise (in channels without strong maser 
features) of ~20 mJy beam" 1 . The synthesized beam was 300 by 275 milliarcseconds (mas) 
at a position angle of —83°. The strong masers in SSV13, some 196 arcseconds to the north- 
northwest, contributed sidelobe emission in some spectral channels. In general, however, 
those channels affected did not have maser emission from IRAS 4 A or 4B. With this angular 
resolution, many spectral features mapped to the same position, although two groups each of 
spatially separated masers were seen toward both IRAS 4A and 4B. The positions of these 
four maser groups, as estimated with Gaussian fits to the strongest maser in each group 
using the AIPS task JMFIT, are listed in Table 1. 

For the continuum observations, after the standard calibration, the narrow band data 
were imaged and used in an iterative self-calibration procedure. The self-calibration solutions 
obtained were applied to the broad band data, which were imaged with natural weighting of 
the u,v data (which increases the sensitivity relative to uniform weighting). The synthesized 
beam for this mode was 360 by 325 mas at a position angle of 77deg. The rms noise in this 
image was ~120 /xJy beam -1 . Unfortunately, the weather was not particularly good, and the 
system temperatures were elevated by about a factor 1.5 - 2.0 over good weather conditions. 
Integrated flux densities at A=1.3 cm for IRAS 4A1 and 4BW were 1.33 mJy and 0.54 mJy, 
respectively. IRAS 4A2 and 4BE were not detected to a 3a peak flux density limit of 0.36 
mJy beam" 1 . 



2.2. VLB A Observations 



We have used the VLBA to obtain four epochs of observation of the water maser emission 
associated with a number of low-mass YSOs. Two subsequent papers w ill present results 



on IRAS 16293-2422 and SVS 13. A previous paper (IClaussen et al.lll998l ) presented results 



as 


Marvel et al. 


(b002), while those of 


Wootten et al. 


( 


1999J) and 


Wootten et al. 



The observations were made with one polarization pair (both right and left circular 
polarization) of baseband channels, each with 4 MHz of bandwidth. This provides 54 km 
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s -1 of total velocity coverage. Correlation was made with 256 spectral channels per baseband 
channel, yielding a velocity width in each spectral channel of 0.21 km s _1 . Although the 
velocity coverage is then only ±27 km s _1 from the syste mic, no published higher velocity 



water masers have been d etected from the IRAS 4 sources (IClaussen et al.lll996l ; iBrand et al. 
20031 : iFuruva et alil2003h . 



Since IRAS 4A and 4B are located close together on the sky (within the primary beam 
of the VLBA antennas), we were able to direct the antennas of the VLBA to a position 
between the two sources and then perform two separate correlation passes with the VLBA 
correlator at the precise positions of the sources without loss of sensitivity. The correlation 
positions coincide with the centimeter continuum positions (see Table 1). Since the data 
were obtained contemporaneously, simply being correlated in separate passes, the two data 
sets share identical time ranges and mid-times of observation. Each of the four observations 
were separated in time by approximately three weeks. This separation was chosen b ased on 
past single dish monitoring observations (IClaussen et al.l Il996l ; IWilking et al.l 1 19941 ) , which 
indicated significant flux changes on this timescale. The dates of the observations, mid-times 
of each epoch and Julian date of the mid-times are given in Table 2 along with the number 
of days between each pair of the observations. 

The data from each epoch were reduced using AIPS. Strong continuum sources were 
observed during each epoch to calibrate both the bandpass response of the receiver systems 
of each antenna as well as allowing determination of any residual delays uncompensated for 
by the correlator model. After residual delay removal, a strong maser feature displaying 
simple structure was selected for determination of delay changes during the observations. 
After the removal of these residual delay rates, the same channel was split from the data set 
for iterative self-calibration. The amplitude and phase self-calibration solutions were then 
applied to all spectral channels and the data were edited for poorly calibrated amplitudes. 

After all calibrations were applied, the data were mapped using the AIPS task IMAGR. 
Spectral cubes were formed with pixel cell sizes of 85 /iarcseconds. The restoring beam 
varied in size, but was typically about 800 /iarcseconds by 400 /iarcseconds with a position 
angle very close to —10.0° east of north. The resolution varied for each epoch depending on 
the details of the self-calibration process and whether or not good self-calibration solutions 
were obtained for given antennas at given times. The exact imaging parameters for each 
epoch and rms noises for typical channels (those channels whose noise floor is not dictated 
by dynamic range constraints) are listed in Table 3. 

To aid in the analysis of the maser motions observed, two-dimensional Gaussian compo- 
nents were fit to each maser spot in each channel where the emission from an individual spot 
exceeded five times the RMS noise for that channel. This was carried out using an automated 
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routine in AIPS (SAD) in an iterative process that first fit the strongest Gaussians in a given 
channel, subtracted the resulting component from the map and then reiterated the search for 
Gaussian components above the noise cutoff in the map. This produces a fairly reliable set 
of Gaussian components without the very tedious process of fitting each component above 
the limit in each channel manually. The resulting positions were then grouped together in 
velocity space and position in tabular form. Only masers detected in two or more channels 
above the cutoff and no more than one beamwidth offset in position were retained. Fewer 
than 10% of the components were excluded in this process, pointing to the reliability of the 
iterative fitting process. 

The surviving component fits were then visually inspected using the 3-dimensional data 
analysis tool XGOBjf). Using XGOBI, it is possible to excise spurious Gaussian fits that result 
from the automated Gaussian fitting routine we utilized. If a fit was called into question by 
the XGOBI inspection process, the original maps were inspected prior to deleting the data 
point. Usually, spurious fits were caused by residual sidelobes being fit in several channels 
by the SAD routine in AIPS and occurred in channels that were dynamic-range limited. The 
components were then used for all subsequent fitting and analysis as described in Appendix 
A. 

3. Description and Analysis of the Maser Regions 
3.1. IRAS 4A 

At the time of our observations, the water masers of IRAS 4A were found in two domi- 
nant regions, separated by 251 milliarcseconds (Epoch I) or 80.3 AU. A line connecting the 
two maser regions lies at a position angle of about —50° (north through east) roughly per- 
pendicular to the CO outflow axis. Comparison with our VLA observations show that the 
masers are associated with IRAS 4A2. Figure 1 shows the spatial and kinematic distribution 
of the masers for all four epochs, each with a closeup of the SE maser component. 

The NW region exhibited only a single maser spot. Its peak flux decreased in time from 
0.97 Jy at the first epoch to 0.12 Jy at the final epoch. The maser is at or slightly red-shifted 
relative to the molecular core velocity of +6.7 km s _1 . The SE region contains more highly 
red-shifted maser spots in Epochs I and II, including the maser used to self-calibrate each 
epoch and another weaker maser just west of the reference spot. In Epochs III and IV, 



3 XGOBI, a multivariate data analysis tool that operates under the XI 1 system is freely available for 
download at http://www.research.att.com/areas/stat/xgobi/ 
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another maser appeared ~10 mas north of the reference feature at a velocity of ~12.5 km 
s -1 (Figure 1). Because we self-calibrated the data using the brightest spot as the reference, 
its position remains constant over time and appears at the origin. 

The separation vector from the reference feature to the NW feature increased in mag- 
nitude over time, expanding from 80.4 AU in the first epoch to 82.3 AU in the final epoch. 
The projected speed of expansion is 53.0 km s -1 given the assumed distance of 320 pc to the 
source. Since the line of sight velocity difference between the two regions is only 2 km s _1 
and assuming that this velocity difference and the proper motion result from the projected 
velocity of an outflow jet, we can calculate an inclination for the outflow of 2°. The true 
space velocity is therefore nearly equal to that measured in the plane of the sky, namely 53 
km s _1 . As shown in Table 4, the proper motion vectors for the NW feature relative to 
the reference feature has a position angle of —49.8°, identical to the position angle of the 
separation vector. 



3.2. IRAS 4B 

Similar to IRAS 4A, IRAS 4B exhibits water maser emission in two dominant well- 
separated regions associated with IRAS 4BW: one to the NW and another to the SE. Masers 
in the NW region are red-shifted relative to the cloud velocity of 7 km s _1 with velocities 
ranging from 13 to 20 km s _1 . Masers in the SE regions are blue-shifted with velocities 
ranging from -3 to 4 km s _1 . The radial velocities for the NW and SE groups are consi stent 



with t hose in the northern and southern lobes of the bipolar HCN outflow observed by I Choi 



(120011 ) and the bipolar H 2 CO outflow mapped by lDi Francesco et al.l (120011 ). A complete list 



of the maser component velocities, fluxes, and positions can be found in Appendix A. 

Figures 2 and 3 show the overall spatial and kinematic distribution of the maser emission 
as observed in each epoch for the two maser groups in IRAS 4B. Over time, the masers in 
the NW expand away from the reference feature at (0,0) in the SE maser group. A line 
connecting the two regions lies at a position angle of roughly —29° and the two regions 
are separated by 488 milliarcseconds (mas) or 156 AU. The SE grouping is roughly arc-like 
with an overall extent of 38 mas (12 AU) and a thickness of 2 mas (0.6 AU), while the 
NW clump is more linear and extends about 26 mas (8.3 AU) with a thickness roughly the 
same as the SE clump. The highest velocity maser emission lies at the northern tip of the 
NW group (+20 km s _1 ) and at the southern tip of the SE arc (—3 km s _1 ). In addition, 
a small maser group is located approximately 100 mas (32 AU) to the NNW of the main 
SE arc. This region's structure changed during the epochs, but usually consisted of two or 
three well-separated maser features (by several mas or about 0.3 AU). This distinct region 
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is located close to the line connecting the NW and SW maser groups. If the active maser 
regions mark the interaction of a bipolar jet from IRAS 4BW with ambient gas, then the 
opening angle of the jet would be ~10°. 

Beginning in the second epoch, but much more clearly shown in the third and fourth 
epochs, a new arc-like ridge of lower velocity emission became visible to the south of the 
other masers in the NW region. This region is not well-characterized by point-like emission, 
but is more diffuse and seems to be an extended ridge of maser emission. By the fourth 
epoch, the entire NW feature had grown to roughly 40 mas (12.8 AU) in extent. 

We compared the positions of the reference feature, which is present in all epochs in the 
SE clump (B-REF in Table 4), with a maser in the NW clump that persisted through all four 
epochs (B-5 in Table 4). The separation vector for the two masers increased over 59.81262 
days from 499.59 mas (159.87 AU) to 508.25 mas (162.64 AU). A linear regression fit to the 
magnitude of the separation vector as a function of time yields an expansion velocity of 78.2 
km s _1 with a formal error of ± 3.9 km s _1 . The radial velocity difference between the two 
masers of 18.5 km s™ 1 suggests an inclination of the outflow of about 13° from the plane of 
the sky. 

To investigate the proper motions of the masers associated with IRAS 4BW, we identified 
all components that were present in all four epochs of observation. Three components were 
identified in the NW and two were identified in the SE. We averaged the mean positions 
of the NW and SE components and averaged these two positions to estimate an average 
geometric fiducial position for each epoch. This position was then subtracted from each 
maser component. These referenced positions are presented along with the fluxes and original 
positions in Table 4. Finally, we fit straight lines to the positions of the components as a 
function of time. The resulting tangential velocities and position angles from these fits are 
also presented in Table 4 and are very similar for all five maser components. Their motions 
are very uniform over the four epochs with average values of 43 ± 2 km s _1 at a position 
angle of —38° ± 2°. Figure 4 shows all of the maser components as observed in Epoch 1 and 
their corresponding proper motion vectors. Again, an overall expansion is apparent along 
the axis connecting the two regions. 

4. Discussion 
4.1. IRAS 4A 

Since the water masers toward IRAS 4A have always been found coincident with the 
position of A2 (when they were detected), we consider three alternatives for the origin of 
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the water maser emission seen using the VLB A: 1) the water masers are in an outflow which 
is the close-in extension of the large-scale molecular outflow, driven by IRAS 4A2; 2) the 
water masers are excited by an interaction of the close-in flow and the accretion disk itself 
and 3) the water maser outflow is driven by a companion to IRAS 4A2, as yet unresolved. 
We consider each of these alternatives in turn. 



4-1.1. Alignment with the large scale outflow 



The alignment of the water maser outflow is curious in that it does not have the same 
position angle as any of the larger-scale outflows emanating from the two sources i n IRAS 4A. 



The axis of the large-scale outflow defined by the SiO, CS, CO, and H 2 emission ( iChoi et al. 



20061 ; IChoill2005l ; iBlake et al.lll995l ). and driven by A2, may be drifting with a rate of 0.011 
degrees/ yr (counterclockwise as seen on the sky) and a zero-point angle of around —166° 
(i.e., presumably the jet-injection angle), depending upon the fit of the large-scale emission 
(see Figure 3a of Choi et al. 2006). For the large-scale flow, the blue-shifted emission is 
generally in the south to southwest and the red-shifted emission is generally in the north 
and northeast. The maser outflow, however, on a scale of 200 mas, has a position angle of 
—50°, and has the more highly red-shifted emission to the southeast. Thus, if the masers 
were marking the close-in extension to the flow, the position angle of the flow would have to 
change from +130° to —165°, in a counterclockwise direction, from 40 to 300 AU from the 
star. We consider that such a shift in jet direction is very unlikely, if not unphysical. 



4-1.2. Interaction with the A2 accretion disk 



The jet direction near th e base of the flow , as indicated by lChoi et al.l (120061 ). and also 
by the CS J=7^6 emission (IBlake et al.lll995l ) appears to be about degrees (i.e. north- 
south). If the masers were excited by an interaction of the outflow with the A2 accretion 
disk itself, then the opening half-angle of the jet at its base would have to encompass the 
direction of the maser axis. Thus the opening half-angle of the jet would be ~ 50 deg. Such 
a large opening angle would seem inconsistent with the Class phase and the high degree 
of collimation seen in the larger scale outflow. 
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4-1.3. An unresolved companion to IRAS 4A2 



Perhaps the simplest explanation for the direction of the axis of the water maser emission 
is that it is not related directly at all to the large-scale molecular outflow, but that the masers 
are located in an independent outflow from another young stellar object. A large fraction 
of the water masers around low-luminosity YSOs, which have been observed and mapped 
with very high angular resolution, show bipolar flows and structural evidence for bow shocks 
as the jet impinges on and sweeps up ambient molecular material. Although the structural 
evidence is sparse (because of the small number of maser spots) for IRAS 4A, the bipolar 
velocity separation st rongly suggests that t he masers are tracing a jet outflow as, for example, 
in IRAS 05413-0104 (IClaussen et al.l 119981 ). It is clear from VLA observations (e.g., Furuya 
et al. 2003; this paper) that the water masers toward IRAS 4A are always found toward 
A2 rather than Al, so Al cannot drive the maser outflow. Therefore we suggest, that, in 
this scenario, there must be a companion to A2 that drives the water masers. Since the 
masers are clearly associated with the position of A2, and not Al, it is clear that the maser 
outflow, in this scenario, is not associated with Al. Indeed, the pr esence of a companion to 
A2 separated by 30-80 AU has been proposed by lChoi et al.l (120061 ) to explain the drifting of 
the 4A2 outflow axis. A companion within 100 mas (32 AU) of A2 would be consistent wi th 
it not being resolved in A=3.6 cm continuum observations with the VLA (jReipurthl 120001 ) . 



4.2. IRAS 4B 

While the masers trace a jet similar in velocity to the HCN outflow, the p osition angle 



of —29° is significantly different from the north-south orientation of the outflow. 'Choil (120011 ) 
suggests a dynamical lifetime for the molecular outflow of ~200 years (corrected to a distance 
of 320 pc). This is much longer than the dynamical lifetime of the maser emission which 
is on the order of a year. We conclude that the jet must be precessing in such a way as 
to move the projected outflow axis from a position angle of 0° to —29° in several hundred 
years. In this picture, the jet would be impacting the walls of the outflow cavity as it 
precessed "clockwise" in the plane of the sky. The NW maser group would represent the 
western wall of the pre-existing cavity and the SE maser group the eastern wall. This 
would explain the linear, arc-like structure of the maser groups and the appearance of new 
masers at the southern extent of the NW group. This model can only be r econc iled with 



observations of dense outflowing gas traced in H 2 CO by IDi Francesco et al.l (120011 ) if there 
were a second outflow from an unresolved companion to IRAS 4BW. If our picture is correct, 
one would expect future VLBI observations to reveal new masers south (north) of the NW 
(SE) group. One would also expect that on a larger scale, tracers of the jet would display 
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a characteristic S-shape, bending back to a north-south orientation. Indeed, s 



jets have been proposed fo r a number of YSO outflows ([Matthews et al.ll2006 



20051 ; iTerquem et al.lll999l ) that suggest the presence of a nearby companion (lTerquemlll998l ). 



ow precessmg 



Hodapp et al. 



app 

3mj[ 



The proper motion vectors for the five persistent maser spots relative to a fiducial point 
(Figure 4) suggests a larger position angle than that of the outflow axis defined by the 
relative motion of SE and NW masers (—38° vs. —29°, see Table 4). This is reminiscent 
of th e proper motions seen in the triple radio continuum source in Serpens (jCuriel et al. 
19931 ). where the outer components of the radio jet have the same tangential velocity but 
are moving in slightly different directions. Although the maser emission in IRAS 4B traces 
the jet at scales 40 times smaller than that traced by the radio continuum in the Serpens 
jet, this difference in relative motions may be suggestive of similar physical processes in the 
two jets. 

The larger position angle defined by the proper motion vectors in IRAS 4B suggest 
that, in addition to the tangential velocity imparted to the masing gas by the outflow jet at 
a position angle of —29°, there is an additional component of about 6 km s" 1 perpendicular 
to the outflow axis. It is possible that this motion is related to the expansion of the outflow 
cavity into the ambient cloud. 



4.3. Maser Geometry 

We estimate that the maser outflows in both IRAS 4A and 4B are nearly in the plane 
of the sky (inclination of only 2° for IRAS 4A and about 13° for IRAS 4B). In Claussen 
et al. (1998) we estimated that the inclination of the maser outflow for IRAS 05413-0104 
was only 4°. Models of maser emission (e.g. Hollenbach 1997) suggest that the masers trace 
shocks produced by the interaction of the jet and ambient material along the jet. Shocks 
that propagate close to the plane of the sky provide much longer maser coherence lengths 
than those that propagate at larger inclination angles. Thus outflows that have inclinations 
close to the plane of the sky provide a more favorable maser geometry, perhaps explaining 
why some well-known outflows from low-mass YSOs do not show maser emission. 



5. Conclusions 

We have observed the water masers associated with IRAS 4A and 4B at VLBI resolutions 
in four epochs over three months. We have determined that the masers are related to the 
jets emanating from these YSOs due to their spatio-kinematic distribution. In both sources 
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the masers are found associated with known single components of multiple systems, with 
total separation velocities between 53 and 78 km s _1 . 

This is further confirmed by the large proper motions measured for both sources, which 
clearly rules out rotation due to the mass constraints placed on the central objects by other 
observations. The water masers of IRAS 4B form arc-like structures roughly 10 AU in length 
and less than 0.6 AU in thickness. The structure of these structures changes rapidly with 
time, with many new maser components appearing and disappearing in just one month. The 
orientation of the structures in the plane of the sky does not agree with larger scale outflow 
angle, which we attribute to a possible unseen very close companion. Future observations 
will have to sample the source more frequently than once every three weeks, with once every 
3-5 days probably providing the best results. 

KBM wishes to thank the AAS Council and AAS Executive Officer, Dr. Robert Milkey 
for allowing him to continue research work in his current position. KBM thanks TEK for 
her continued support. And special thanks to T. Beasley for exposing KBM to the rigors 
of Outback observing and a grueling exercise regimen, during which time a portion of this 
paper was written. Finally, KBM thanks the NRAO for hosting him during a research leave 
from the AAS, during which time the bulk of this paper was completed. 

A. Table of Maser Components for IRAS 4B 

Maser component fits were averaged using a flux-squared weighting scheme, that em- 
phasizes the strongest portion of the maser emission and produces a reliable location for the 
maser. As some masers are present for only three channels, a full Gaussian fit to the maser 
as a function of velocity was deemed impractical, although it could have been used for the 
masers that existed over many channels. We opted to use the same position-determining 
scheme for all features. These averaged component fits and their corresponding errors are 
given in Table 5 for IRAS 4B. Note that the uncertainty in position for the reference channel 
is simply the error on the mean position values determined for the spot and is indicative of 
the reliability of the Gaussian fitting routines in the presence of noise in the images. 
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Table 1: IRAS 4 A and IRAS 4B Properties 



Property 






Value 




Reference 


IRAS4A 












Al 2.7 mm Position (J2000.0) 


3 h 29m 


10 


s .510 +31° 13'31" 


.01 


Looney, Mundy and Welch 2000 


Al 1.3 cm Position (J2000.0) 


3 h 2 gm 


10 


s .527 +31° 13'31" 


.06 


This study 


Al 3.6 cm Position (J2000.0) 


3 h 2 gm 


10 


s .529 +31° 13'31" 


.05 


Reipurth et al. 2002 


A2 2.7 mm Position (J2000.0) 


3 h 29 m 


10 


s .413 +31° 13'32" 


.20 


Looney, Mundy and Welch 2000 


A2 3.6 cm Position (J2000.0) 


3 h 2 gm 


10 


s .421 +31° 13'32" 


.21 


Reipurth ct al. 2002 


A2 H 2 NW Position (J2000.0) 


3 h 2 gm 


10 s 


.4101 +31° 13'32" 


.235 


VLA-12Aug98; This study 


A2 H 2 SE Position (J2000.0) 


3 h 2 gm 


10 s 


.4226 +31° 13'32" 


.142 


VLA-12Aug98; This study 


Al Peak Flux (2.7 mm) 






107 mjy 




(2") Looney, Mundy, & Welch 2000 


A2 Peak Flux (2.7 mm) 






23 mjy 




(2") Looney, Mundy, & Welch 2000 


Al Peak Flux (3.6 cm) 






0.32 mjy 




Reipurth et al. 2002 


A2 Peak Flux (3.6 cm) 






0.11 mjy 




Reipurth et al. 2002 


Bolometric Luminosity (Lq) 






14 




Andre, Ward-Thompson, & Barsony 2000 


Bolometric Temperature (K) 






34 




Andre, Ward-Thompson, & Barsony 2000 


IRAS4B 












2.7 mm Position (J2000.0) 


3 h 2 gm 


11 


s .988 +31° 13' 08" 


.10 


Looney, Mundy, & Welch 2000 


1.3 cm Position (J2000.0) 


3 h 2 gm 


12 


s .001 +31° 13' 08" 


.16 


This study 


3.6 cm Position (J2000.0) 


3 h 2 gm 


12 


s .003 +31° 13'08" 


.14 


Reipurth et al. 2002 


B H 2 NW Position (J2000.0) 


3 h 2 gm 


11 s 


.9923 +31° 13'08" 


.363 


VLA-12Aug98; This study 


B H 2 SE Position (J2000.0) 


3 h 29m 


12 s 


.0104 +31° 13'07" 


.955 


VLA-12Aug98; This study 


Peak Flux (2.7 mm) 






58 mjy 




(2") Looney, Mundy, & Welch 2000 


Peak Flux (3.6 cm) 






0.33 mjy 




Reipurth ct al. 2002 


Bolometric Luminosity (Lq) 






14 




Andre, Ward-Thompson, & Barsony 2000 


Bolometric Temperature (K) 






36 




Andre, Ward-Thompson, & Barsony 2000 
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Table 2. Dates and times of Observations 



Epoch 


Date 


Mid-time 


Julian Date 


ATime [d] 




of 1998 


UT 


-2451000.0 




I 


September 15 


11:01:23 


71.95928 


0.00000 


II 


October 4 


10:01:30 


90.91771 


18.95843 


III 


October 27 


08:01:30 


113.83437 


41.87509 


IV 


November 14 


06:31:30 


131.77187 


59.81259 
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Table 3. Imaging parameters for each epoch 



Epoch 


Beam Dimensions 


Beam Position Angle 


RMS 




[//arcsec] 


[° E of N] 


[mJy/Beam] 


4A-I 


1260 x 910 


11.8 


10.4 


4A-II 


830 x 440 


-21.2 


6.8 


4A-III 


810 x 400 


-19.1 


11.7 


4A-IV 


650 x 370 


-8.2 


9.6 


4B-I 


650 x 350 


-14.4 


13.9 


4B-II 


620 x 330 


-12.4 


7.6 


4B-III 


590 x 330 


-10.8 


9.9 


4B-IV 


720 x 420 


-5.8 


6.0 
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Table 4. Proper Motions for Masers Observed in All Four Epochs 



Epoch ID Velocity 


Peak Flux X-Offset 


X re / 


Y-Offset 




(km sec -1 ) 


(Jy) (mas) 


(mas) 


(mas) 


(mas) 



IRAS 4A 



I-A-REF a 


10.40 


3.29 


0. 


.083 






0.073 






II-A-REF 


10.55 


4.14 


-0. 


063 






-0.185 






TTT A RF,F 


1U.UO 


7 Q7 


n 

u . 


DD7 






u.uou 






IV-A-REF 


10.72 


14.9 


-0. 


005 






0.002 






I-A-l 


7.61 


0.969 


-167. 


.543 


-167. 


.626 


187.223 


187. 


.150 


II- A- 1 


7.65 


0.211 


-169 


.013 


-168. 


.950 


187.985 


188. 


.170 


III-A-1 


7.71 


0.208 


-170. 


,631 


-170. 


.638 


189.606 


189. 


,636 


TV A 1 

1 V jV 1 


7 8^ 


n 19^ 


1 79 


.uyo 


1 79 


,uyo 


1 QD Q94 


1 Qfl 


Q91 


Vexp 


53.02 


















PA(Deg) 


-49.82° 


















I-A-2 


10.34 


1.90 


-1. 


.859 


-1. 


,942 


-0.086 


-0 


.160 


II-A-2 


10.36 


2.26 


-1. 


906 


-1. 


,844 


-0.118 


0. 


,067 


III-A-2 


10.29 


1.74 


-1. 


.844 


-1. 


,851 


0.695 


0. 


,725 


IV-A-2 


10.35 


1.27 


-1. 


.938 


-1. 


933 


1.113 


1. 


,108 


IRAS 4B 




















I-B-REF b 


-0.750 


36.8 


-0. 


.020 


122. 


,668 


0.0782 


-227 


,845 


II-B-REF 


-1.161 


20.4 


-0. 


.001 


123. 


,598 


-0.041 


-228. 


,993 


III-B-REF 


-0.224 


56.6 


-0. 


.000 


124. 


,475 


-0.013 


-229 


,284 


IV-B-REF 


-0.273 


34.3 


-0. 


002 


125. 


,498 


0.010 


-230 


.661 


I-B-l 


0.380 


0.368 


-1. 


736 


120. 


932 


10.016 


-217. 


.830 


II-B-1 


0.302 


2.14 


-1. 


.677 


121. 


,921 


9.384 


-219 


,549 


III-B-1 


0.485 


0.256 


-1. 


.626 


122. 


,849 


8.726 


-220 


.559 


IV-B-1 


-0.158 


0.09 


-1. 


732 


123, 


,767 


8.912 


-221 


,749 



PA(Deg) 144.12' 
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Table 4 — Continued 



Epoch ID Velocity Peak Flux X-Offset X re/ Y-Offset Y re/ 
(km sec -1 ) (Jy) (mas) (mas) (mas) (mas) 



I- B-2 2.81 0.829 

II- B-2 3.05 0.134 

III- B-2 3.25 1.24 

IV- B-2 3.39 0.218 
v exp 41.03 

PA(Deg) 140.07° 

I- B-3 17.46 0.634 

II- B-3 17.07 2.80 

III- B-3 17.27 0.521 

IV- B-3 17.13 0.318 
Vexp 41.90 

PA(Deg) -40.50° 

I- B-4 17.74 2.31 

II- B-4 17.66 5.90 

III- B-4 17.60 1.50 

IV- B-4 17.43 0.862 
V e xp 41.47 

PA(Deg) -36.51° 

I- B-5 17.95 0.643 

II- B-5 17.99 0.957 

III- B-5 17.84 1.24 

IV- B-5 17.22 1.55 
Vexp 46.15 

PA(Deg) -36.44° 



-9 


.429 


113 


.239 


33 


.556 


-194. 


290 


-9 


.531 


114 


,067 


33 


.844 


-195. 


.089 


-9 


.074 


115. 


.402 


32 


,605 


-196. 


.679 


-9 


.417 


116 


.082 


32 


,975 


-197. 


.686 


-239 


.832 


-117 


,165 


432 


,737 


204. 


.892 


-241 


.622 


-118. 


.024 


434. 


.953 


206. 


020 


-243 


.675 


-119. 


.200 


436. 


.448 


207. 


.164 


-245. 


.662 


-120 


.163 


439. 


.063 


208. 


.402 


-239 


.610 


-116 


,943 


430. 


,618 


202. 


.772 


-241. 


.452 


-117. 


.853 


432 


,924 


203. 


.991 


-243 


.359 


-118. 


.884 


434. 


,451 


205. 


.167 


-245. 


.111 


-119. 


.612 


437. 


,038 


206. 


.377 


-239 


.816 


-117. 


.148 


438. 


.361 


210. 


.516 


-241. 


.704 


-118. 


.106 


440. 


.878 


211. 


.947 


-243 


.767 


-119 


.292 


442. 


,811 


213. 


.527 


-245. 


.495 


-119. 


.997 


445. 


,035 


214, 


.374 
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a Strongest maser in IRAS 4A used as both a position and velocity reference. 

b Strongest maser in IRAS 4B used as a position reference only. Proper motions 
are measured relative to a fiducial point for each epoch given as X re j and Y re j for 
this maser (see text). 
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Table 5. Fitted Gaussian Component in IRAS 4B 



Epoch ID 


Other 


Velocity 


Peak Flux 


X- Offset 


Y- Offset 


Ox 


Oy 






(km s' 1 ) 


(Jy) 


(mas) 


(mas) 


(/ias) 


(fxa,s) 


T QTT 1 


R 
rS-Z 


o si 
Z.ol 


U.oo 


O /I98 

-y .4Zo 


qq ccc 
00.000 


A A 
4.4 


8 O 
O.U 


T Q"FT 9 




o qo 
Z.OU 


8 9/1 
O.Z4 


7 IPO 

- ( .ooz 


97 8Q1 
Z ( .sy 1 


1 
U. 1 


1. 1 


T C!"R 1 
L-oHi-O 




9 91 
Z.Z1 


4.y i 


K 799 
-O. ( zz 


91 8^8 


1 1 
1. 1 


1 1 
1.0 


T Q"FT A 




9 97 
Z.Z ( 


n 7P. 

U. ( 


K QPP 

-o.yoo 


99 <A9K 


K K 
0.0 


K 1 
0.0 


T Q"FT K 




9 i n 
Z. 1U 


n aa 

U.44 


P A 11 
-0.400 


01 /1/18 
Z0.44O 


7 8 
( .0 


19 

iz.y 


T Q"FT P 




1 Q/1 

i.y4 


n 7P. 

U. ( 


K 107 
-0.0Z ( 


90 Q08 

zu.yuo 


Q 8 
O.o 


K K 
0.0 


T QTT 7 




1 88 
l.oo 


1 7P 
1.(0 


K 88*? 
-O.ooO 


99 m p 

ZZ. UIO 


1 9 
l.Z 


9 K 
Z.O 


T Q"FT 8 




1 7/I 
1.(4 


n c;8 

U.Oo 


K 9/17 
-0.Z4 ( 


90 7^0 

zu. ( oy 


1 

0. 1 


118 
11. o 


T Q"FT Q 




1 .oo 


n 17 
U.O ( 


A Kl K 
-4.010 


1 8 A 71 
lo.4 ( 


7 

( .0 


TAP 
14.0 


T Q"FT 1 n 
1-OI1/-1U 




1 9Q 

i.zy 


U.oo 


Q 1 0/1 
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Fig. 1.— The spatial and kinematic distribution of the water masers associated with 
IRAS 4A2 for the four epochs observed. The top frame in each panel shows the veloc- 
ity and flux of the detected features, the lower frame shows the spatial distribution. The 
maser to the N W is colored blue (at the cloud velocity) and those to the SE are coded red 
(red-shifted relative to the cloud velocity). The size of the circles in the figure is propor- 
tional to the flux of that maser component. All components measured in each channel are 
presented (e.g., no component averaging has been performed). The error in the position of 
the spots (typically less than 20 ^uarcseconds) cannot be plotted on the large scale figure. 
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Fig. 2.— Spatial and kinematic distribution of the blueshifted masers southeast of 
IRAS 4BW for all four epochs. In each panel, the top frame shows the velocity (color- 
coded) and the flux (symbol size) of the detected features, and the lower frame shows the 
spatial distribution and the velocity distribution as a function of declination. The uncer- 
tainty in the fitted position of each feature circle cannot be displayed at this scale, but in 
the worst case is of order 20 /zarcseconds. 
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Fig. 3.— Spatial and kinematic distribution of the redshifted masers northwest of 
IRAS 4BW for all four epochs. In each panel, the top frame shows the velocity (color- 
coded) and the flux (symbol size) of the detected features, and the lower frame shows the 
spatial distribution and the velocity distribution as a function of declination. The uncer- 
tainty in the fitted position of each feature circle cannot be displayed at this scale, but in 
the worst case is of order 20 /zarcseconds. 
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Fig. 4. — A plot of the proper motions for the five maser components present in all four 
epochs for IRAS 4BW relative to their positions in Epoch 1 (see Table 4). The length of the 
proper motion vectors are 2 km s _1 per plot unit. Note that the proper motions are not at 
the same position angle as a line connecting the two maser groups (see text). 



